Error tolerance of the BosonSampling model for linear optics quantum computing 
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Linear optics quantum computing (LOQC) is a promising approach to implementing scalable 
quantum computation (QC). However, this approach has very demanding physical resource require- 
ments. Recently, Aaronson & Arkhipov showed that a simplified model, which avoids the requirement 
for fast feed-forward and post-selection, while likely not capable of solving BQP-complete problems 
efficiently, can solve an interesting sampling problem, believed to be classically hard. Loss and 
mode-mismatch are the dominant sources of error in such systems. We provide evidence that even 
lossy systems, or systems with mode-mismatch, are likely to be classically hard to simulate. This 
is of practical interest to experimentalists wishing to demonstrate such systems, since it suggests 
that even with errors in their implementation, they are likely implementing an algorithm which is 
classically hard to simulate. Our results also equivalently apply to the multi-walker quantum walk 
model. 



Introduction — Quantum computing (QC) !lj offers 
exponential increases in computing power for certain key 
problems, however the building of a full-scale universal 
quantum computer remains well beyond current tech- 
nological abilities. Nevertheless, single purpose machines 
capable of solving particular problems, may become pos- 
sible on a shorter timescale. For example, the problem of 
building a quantum emulator - a well controlled quantum 
system whose dynamics approximate those of a classically 
intractable physical quantum system of interest - may be 
solvable in the medium term. 

Linear optics QC (LOQC) [5J|3] has emerged as a lead- 
ing candidate for the implementation of QC. In princi- 
ple, universal QC operations can be implemented in sev- 
eral different ways using linear optics, photon production 
and counting, quantum memory and fast feedforward 
[4j. Recently Aaronson & Arkhipov (AA) [5] studied an 
alternate, but more straightforward approach to imple- 
menting LOQC. Here only photon preparation, passive 
linear optics elements and photo-detection are required. 
No fast feed-forward or quantum memory is necessary. 
AA showed that while this model probably cannot solve 
BQP-complete problems (i.e. it is likely not universal for 
QC), it can solve a sampling problem closely related to 
the #P-complete matrix permanent problem, believed 
to be classically hard. That is, a passive linear optical 
network, fed with single photon states, can emulate a 
classically hard Boson-sampling problem. 

An unanswered problem with quantum emulation is 
the effect of errors. In a universal QC, errors can be 
corrected via the implementation of fault tolerant codes. 
However, typically the restricted set of operations avail- 
able to an emulator may be insufficient for fault tolerant 
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error correction. Thus a pressing question for particular 
emulator proposals is whether they can still simulate clas- 
sically hard dynamics in the presence of realistic levels of 
noise. 

In this paper we address this problem with respect to 
the AA scheme and its tolerance to loss, i.e. detector, 
source and component inefficiency. This is the major er- 
ror source for optical quantum processing. We present 
strong evidence that the AA scheme can still simulate 
classically hard problems in the presence of realistic inef- 
ficiencies. The trade-off compared to an error-free system 
is an increase in circuit size and photon number. 

The Model — The AA model solves the so-called 
BosonSampling problem by sampling the photon num- 
ber configuration statistics at the output of a linear op- 
tics circuit with a multi-photon input state. We begin by 
preparing some number of modes N, where some config- 
uration of n < N modes are initialised with the single 
photon state |1), and the remaining N — n modes are in 
the vacuum state |0). Thus the input state is of the form 
|^i„) = |li, ... , l n ,0„+i, . . . , CM, where N = 0{n 2 ). The 
input state then passes through a passive linear op- 
tics network consisting only of beamsplitters and phase 
shifters, implementing a unitary map U on the photon 
creation operators, a\ — > Ylj UijOj, where a\ is the cre- 
ation operator acting on mode i. Importantly, Reck et al. 
[6] showed that an arbitrary U can always be decomposed 
into a polynomial number of optical elements. Thus any 
U of this form can always be efficiently experimentally 
constructed. 

In the occupation number representation the output 
state is of the form 

\i> out )=Y / ls\n{ S \ni S) ...,nf) (1) 
s 

where S are the different configurations of photons across 
the N modes in the output state and 7s are the corre- 
sponding amplitudes. The model is illustrated in Fig. [T] 
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FIG. 1: The model for the BosonSampling problem. We be- 
gin with n photons, one in each of some number of modes, and 
none in the remaining modes. After the linear optics network, 
U, we measure some configuration of photons S. 

By detecting photon number at the output the Boson 
distribution for some particular input configuration and 
unitary scatterer, U, can be sampled. AA propose that 
if U is picked randomly then it is not possible for a clas- 
sical computer to efficiently simulate the sampling of the 
distribution. 

We can formulate the problem in terms of a competi- 
tion between two participants, Alice and Bob. Alice pos- 
sesses a quantum emulator whilst Bob possesses a classi- 
cal simulator. The adjudicator, Victor, draws a random 
unitary, U, from a hat and Alice and Bob race to sam- 
ple its distribution to some predetermined accuracy in 
the shortest possible time. According to AA's thesis, for 
sufficiently large systems, Alice will consistently win the 
race. 

Why Boson sampling is hard — AA demonstrate that 
should BosonSampling be classically easy to simulate, 
this would have very surprising consequences in compu- 
tational complexity theory, and therefore suggest that 
it is most probable that BosonSampling is not clas- 
sically easy to simulate [5]- From the complexity-class 
point of view an intuitive reason suggesting that Boson- 
Sampling is classically hard is that each 75 is a func- 
tion of a matrix permanent, P(S) oc Per(Ag). The matrix 
permanent problem resides in the class #P-complctc, 
strongly believed to be classically hard. Importantly 
though, BosonSampling does not allow us to calculate 
matrix permanents efficiently, since doing so would re- 
quire an exponential number of measurements. Rather, 
it samples across many different matrix permanents. 

More pragmatically, we can look at particular classical 
simulation approaches and see how they fail. Considering 
the size of the Hilbert space it can easily be seen that, 
with indistinguishable photons, the number of parame- 
ters \S\ = ('' V+ ™~ 1 ) grows exponentially with the number 
of photons and modes in the system. The emphasis on in- 
distinguishable photons is important. If the photons are 
distinguishable then the output can be shown to factor 
into: 

l«=nE^hif ) 5 4f ) ... 1 n^' ) ) (2) 

i S' 

where now 70/ are the amplitudes for the different out- 
puts a single input photon can reach, with the product 



over all n input photons. For each photon there are only 
S' = N output configurations so the total number of 
independent parameters is n x N. As this only grows 
polynomially for distinguishable input photons, such a 
situation is classically simulatable. We will return briefly 
to this issue later when we consider mode-matching. AA 
suggest that a 400 mode interferometer fed with 20 sin- 
gle photons lies around the boundary of the simulation 
powers of current classical computers. For this case we 
find \S\ « 7 x 10 33 . 

The problem might also be simulated in the quadra- 
ture field basis. It is known that such a simulation will 
be efficient if the input state is Gaussian |5] and the 
photon counting is not adaptive 1 . Single photon states 
are non-Gaussian but have a finite overlap with Gaus- 
sian states. This motivates the question: how close can 
a Gaussian simulation be to the BosonSampling prob- 
lem? In particular, suppose we replace each single photon 
state at the input with a Gaussian state and then simu- 
late the unitary evolution of interest with this state. How 
similar can the two output states be? 

A measure of the similarity of two density operators is 
the trace distance, D = |tr|p — a\. The trace distance 
is a metric. It equals zero when the states are identi- 
cal and one when they are orthogonal. It is monotoni- 
cally related to the lower bound of the Shannon mutual 
distinguishability between the two density operators [3]. 
Thus, given the symmetry of the physical situation we 
consider 2 , the trace distance can bound the amount of 
information that can be learned about density operator 
p via a simulation which produces density operator a. 
Furthermore, the trace distance is invariant under uni- 
tary transformations, thus by computing the minimum 
trace distance between the input single photon distri- 
bution and a Gaussian input state we can evaluate the 
similarity between the quantum emulation and a classical 
simulation of the output for arbitrary unitaries. 

The trace distance between a single photon state and 
single Gaussian state is easily calculated and is minimized 
by maximizing the absolute value of the single photon 
amplitude of the Gaussian. We obtain D — 0.522 for the 
squeezed state with amplitude /3 = v2 and squeezing of 
V = 1/3 (where the quantum noise limit is V = 1) |10) . 
For this simple case of pure states the trace distance be- 
tween n copies of the single photon state and n copies 
of the minimizing Gaussian state is D n = 1 — (1 — D) n , 
which asymptotes to unity exponentially in the number 



1 An adaptive measurement scheme is one in which the type of 
measurement made on the (i+l)th mode is a function of the mea- 
surement result obtained from the ith mode. In a non-adaptive 
measurement scheme, such as in the AA proposal, all the mea- 
surement bases are predetermined. 

2 The density operators are symmetric in the sense that we assume 
every single photon state input in p is replaced by the optimal 
Gaussian approximation in a. Thus the distortion of a away from 
p is symmetric over all the modes. 
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of copies. Thus we find that a Gaussian simulation of 
the problem will rapidly become useless as the size of 
the system grows. Considering again the case of 20 in- 
put photons, and using the above distance minimising 
classical state, we find D n = 1.0 - 3.86 x 10" 7 - clearly 
very close to one, implying that the Gaussian simulation 
would give virtually no information about the BOSON- 
Sampling problem. 

Effect of Loss — Having gained some physical intu- 
ition into why the BosonSampling problem is hard we 
now wish to investigate the effect of errors, in particular 
photon loss, on the operation of the quantum emulation. 
In Appendix A we present an argument based on trunca- 
tion of the photon number basis in a classical simulation 
to show the problem becomes "easier" as loss increases. 
Nevertheless, if only small amounts of loss are present 
then it may still be reasonable to post-select successful 
events in a probabilistic way. Provided the overhead is 
not too large the AA system can still be said to be solv- 
ing a hard problem efficiently. 

Consider the case of the 20 input photons to the 400 
mode interferometer. Let us assume the total efficiency 
for a single photon to be produced, transmitted through 
the interferometer and to be successfully detected is rj. 
Let us further assume that this efficiency is the same 
for all photons. Under such conditions the detection and 
transmission efficiencies can be commuted back to the 
source and lumped together as a single source efficiency 
[TT1 [12] . The chance that all 20 photons will be detected 
at the output is then P — r/ 20 . Provided this number 
is not too low then statistics can still be obtained for 
the sampling problem by running the system many times 
and only keeping those cases where 20 photons are actu- 
ally detected. Suppose we conservatively require P = 0.5, 
meaning we need to run the system about twice as many 
times as the no loss case to obtain good statistics. We 
conclude that a bulk efficiency of about i] ^ 0.96 is 
needed. Lower efficiencies are allowed if we allow for more 
repetitions of the experiment, however if Alice postselects 
too much she will lose her advantage over Bob's classical 
processor. The figure of 0.96 is good compared to thresh- 
olds for universal quantum computing where efficiencies 
higher than 0.99 are typical, and clearly illustrates the 
advantage of single purpose machines. However, this still 
represents a major challenge given current technology. 
In the following we consider a different direction, more 
promising for the short term. 

Is lossy Boson sampling a hard problem? — We now 
address the question 'can lossy AA systems be efficiently 
classically simulated?', and present evidence that this is 
not the case in situations of near term experimental in- 
terest. Although we have seen that if the efficiency is too 
low then the probability of post-selecting the samples of 
a lossless scattering problem are very low, it may still be 
true, for moderate levels of loss, that the lossy Boson- 
Sampling problem itself 'is hard to simulate on a classical 
computer - and thus of interest for quantum emulation. 
We base our evaluation on whether the problem remains 



hard for the two simulation approaches discussed earlier. 
We observe that if the average photon number of the 
input state is equal to, or greater than that for the loss- 
less case, with an equal total number of modes, then the 
number of parameters, \S\, required for a direct number 
state simulation will be at least as large as for the loss- 
less case. On the other hand we recognize that a lossy 
single photon state is more Gaussian than a pure one, 
and so we could expect a Gaussian simulation of lossy 
BosonSampling to be a better approximation than it 
would be for lossless BosonSampling. However, we also 
know from the previous section, that for small amounts of 
loss the AA system remains hard to simulate due to the 
ability to post-select. By calculating the trace distance 
between the output state of an AA system with a small 
amount of loss and the corresponding closest Gaussian 
state we can obtain a benchmark condition for deciding 
if a Gaussian simulation will be useful in more general 
situations. In summary we seek lossy A A systems that 
fulfil the following criteria: 

(i) nL ^ n, where t%l (n) is the average photon number 
for the lossy (lossless) case and the total number of modes 
is of N = 0{n 2 ); and 

(ii) Dl ^ D, where Dl is the trace distance between 
the output of a lossy AA system and the closest equiv- 
alent Gaussian state, whilst D is similar but for an AA 
system for which the loss is sufficiently small that the 
post-selection probability P ^ 0.5. We particularly con- 
sider the case of n = 20 for our calculations. 

A lossy single photon state can be represented as a den- 
sity operator of the form pi ossy = (1 — 7?)|0)(0| + 
In general the closest Gaussian state will be a phase dif- 
fused, displaced squeezed state. The expression for the 
trace distance between n lossy single photon inputs and 
n phase diffused, displaced squeezed state becomes: 
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is the ith photon number probability of a phase diffused, 
displaced squeezed state with displacement f3 and squeez- 
ing V = {ijl — v) 2 (with /i 2 — v 2 = 1) , and Hi are Hermite 
polynomials. The trace distance can be minimized as a 
function of f3 and V to give -D™ ln . In particular we find 
that for n = 20 and r\ = 0.96 that D™ in = 0.9987 for 
= 1.217 and V = 0.3801. 

In Fig. [2jtop) we plot -D™ ln against the loss rate rj for 
various input photon numbers, n. Two trends are clear. 
For a fixed input photon number, increasing loss leads 
to decreasing minimum trace distance, indicating that 
the Gaussian simulation is becoming a better approxi- 
mation to the quantum emulation. On the other hand, 
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for fixed loss, increasing the input photon number in- 
creases the minimum trace distance, indicating that the 
Gaussian simulation is becoming a worse approximation 
to the quantum emulation. This suggests that the "hard- 
ness" of the simulation can be maintained in the presence 
of increasing loss by simultaneously increasing the input 
photon number. 

In Fig. [^middle) we plot the number of photons that 
must be input into the system such that the average pho- 
ton number after loss, t%l — m, and hence the number 
of configuration amplitudes that would need to be cal- 
culated via a number state simulation remains approx- 
imately constant. Again we find that the "hardness" of 
the simulation can be maintained if we compensate for 
loss by injecting additional photons. 

In Fig. [^bottom) we plot the regions satisfying our 
conditions (i) and (ii), i.e. the regions where the "hard- 
ness" of the simulation, with respect to photon number 
and Gaussian approaches, is maintained at the same level 
as for n = 20, rj = 0.96. We find that for large pho- 
ton numbers the requirements for a hard calculation are 
rather lenient. In particular, loss levels can exceed 50% 
for n > 200. 

Mode-mismatch — In addition to loss, mode-mismatch 
is another dominant source of error in AA type systems. 
Using the techniques of Ref . [T3] , we can decompose the 
spectral distribution function of the input photons into 
overlapping and non-overlapping components, governed 
by a distinguishability parameter A. As A grows the num- 
ber of independent configurations drops as the number of 
distinguishable photons increases. 

The standard technique to overcoming mode-mismatch 
is using narrowband filtering. This technique is widely 
used in virtually all LOQC experiments. The idea here is 
that if we project onto a narrow window in say frequency- 
space, then within that window distinguishable photons 
will appear indistinguishable. The introduction of filter- 
ing is equivalent to loss, since we are discarding a por- 
tion of the wave-function. Specifically, after filtering onto 
some small window w around uiq, the probability of de- 
tecting the photon is p(wo) — f"° W \ip(u))\ 2 dui, where 
ip{uj) is the spectral distribution function of the photons. 
This is equivalent to a lossy channel with loss probabil- 
ity 1— p(wo). Thus, the problem of mode-mismatch is 
reduced to the problem of loss. Our results then indicate 
that some level of mode mismatch can be tolerated pro- 
vided the total loss budget, including filtering, does not 
exceed the amounts allowed by Fig. [^bottom). 

Relation to multi-walker quantum walks — Quantum 
walks (QWs) [TIHTT] have recently obtained a lot of in- 
terest as an approach to quantum information process- 
ing. It was noted in Ref. [18] that to achieve exponential 
speedup in quantum algorithms, multiple walkers (pho- 
tons) must be introduced into a QW system (assuming 
non-exponentially sized graphs). A formalism for multi- 
walker walks was introduced by Rohde et al. |19j . In- 
deed multi-walker systems have begun to be experimen- 
tally demonstrated [501 121] • The difference between gen- 
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FIG. 2: Colour online, (top) Minimum trace distance between 
the density operators of the actual lossy system compared to 
a Gaussian simulation using an optimised squeezed state for 
different input photon numbers n. (middle) Physical resource 
requirements, as measured by the required number of photons, 
such that the average number of photons in the system is m. 
(bottom) The region satisfying our conditions (i), in blue, and 
(ii), in red. The overlapping region is presumed to be hard to 
classically simulate. 



eral linear optics networks and the QW formalism is the 
presence of coins in quantum walks. However, it can be 
shown |22j that a multi-walker quantum walk can be di- 
rectly mapped to a beamsplitter network with quadratic 
resource overhead. A reverse mapping also always exists. 
Thus, QWs can be regarded as isomorphic to the formal- 
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ism presented by AA, and the results we present here 
apply equally to the multi-walker QW formalism (in the 
absence of feed- forward) . Furthermore, AA suggest that 
it is unlikely that general LO networks are universal for 
QC, implying the multi- walker QW formalism is also not 
universal. 

Conclusion — We have considered the emulation of 
the BosonSampling model with LOQC and investi- 
gated the effect of loss and mode-mismatch. We have 
presented evidence that lossy systems remain hard to 
simulate provided the loss is compensated through the 
injection of higher photon number input states. This ob- 
servation should be of interest to quantum optics exper- 
imentalists, for whom experimental implementations of 
AA systems are in principle viable, but where loss is a 



ubiquitous problem. These results suggest that even with 
significant loss in their systems, they may be able to ex- 
perimentally implement a quantum emulation in optics 
which could consistently beat the best classical simula- 
tion. 



Acknowledgments 

This research was conducted by the Australian Re- 
search Council Centre of Excellence for Quantum Com- 
putation and Communication Technology (Project num- 
ber CE110001027). 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

[2] E. Knill, R. Laflamme, and G. Milburn, Nature 409, 46 
(2001). 

[3] P. Kok and B. W. Lovett, Introduction to Optical Quan- 
tum Information Processing (Cambridge Press, 2010). 

[4] P. Kok, W. J. Munro, K. Nemoto, T. C. Ralph, J. P. 
Dowling, and G. J. Milburn, Rev. Mod. Phys. 79, 135 
(2007). 

[5] S. Aaronson and A. Arkhipov, STOC (2011), 

arXiv:1011.3245. 
[6] M. Reck, A. Zeilinger, H. J. Bernstein, and P. Bertani, 

Phys. Rev. Lett. 73, 58 (1994). 
[7] S. D. Bartlett and B. C. Sanders, Phys. Rev. Lett. 20, 

207903 (2002). 

[8] S. D. Bartlett and B. C. Sanders, Phys. Rev. Lett. 88, 
097904 (2002). 

[9] C. A. Fuchs and J. van de Graaf (1997), quant- 

ph/9712042. 
[10] H. P. Yuen, Phys. Rev. A 13, 2226 (1976). 
[11] A. P. Lund, T. C. Ralph, and H. L. Haselgrove, Phys. 

Rev. Lett. 100, 030503 (2008). 
[12] M. Varnava, D. E. Browne, and T. Rudolph (2007). 
[13] P. P. Rohde, W. Mauerer, and C. Silberhorn, New J. 

Phys. 9, 91 (2007). 
[14] Y. Aharonov, L. Davidovich, and N. Zagury, Phys. Rev. 

A 48, 1687 (1993). 
[15] D. Aharonov, A. Ambainis, J. Kempe, and U. Vazirani, 

in Proceedings of the 33th STOC, New York 50 (2001). 
[16] A. M. Childs, Phys. Rev. Lett. 102, 180501 (2009). 
[17] J. Kempe, Cont. Phys. 44, 307 (2003). 
[18] P. P. Rohde (2010), arXiv: 1010.4608. 
[19] P. P. Rohde, A. Schreiber, M. Stefanak, I. Jex, and C. Sil- 
berhorn, New J. Phys. 13, 013001 (2011). 
[20] A. Peruzzo, M. Lobino, J. C. F. Matthews, N. Matsuda, 

A. Politi, K. Poulios, X.-Q. Zhou, Y. Lahini, N. Ismail, 

K. WorhofT, et al., Science 329, 1500 (2010). 
[21] J. O. Owens, M. A. Broome, D. N. Biggerstaff, M. E. 

Goggin, A. Fedrizzi, T. Linjordet, M. Ams, G. D. 

Marshall, J. Twamley, M. J. Withford, et al. (2011), 

arXiv:1103.0604. 



[22] P. P. Rohde, A. Fedrizzi, and T. C. Ralph (2011), 
arXiv: 1102.4385. 



I. APPENDIX 

A. Simulation via truncation 

We consider an alternate approach to classically simu- 
lating a lossy AA system: via truncation of the classical 
simulation to a given number of photons and discarding 
higher photon number terms. Defining an error metric as 
the distance between the probability distributions for the 
ideal and truncated simulations, we find that 

6 = E -PlossY [1 - (1 - 1/NY] , (5) 

i=m+l ^ ' 

where m is the number of photons the simulation is trun- 
cated to, and we have made the assumption that the net- 
work is approximately balanced (i.e. an incident photon 
has approximately the same probability of reaching any 
given output). 

Fig. |2jbottom) illustrates the error measure against 
the number of photons the simulation is truncated to. 
When there is no loss in the quantum system, the clas- 
sical simulation is a very poor approximation unless we 
simulate the full n photons, which is presumed to be ex- 
ponentially hard via the results of AA. As the loss rate 
increases we can tolerate greater levels of truncation in 
the classical simulation, making it more tractable. 

As before, when the loss rate is very high, we are ef- 
fectively simulating the vacuum state, which is compu- 
tationally easy. However, for moderate loss rates, where 
there are still a large number of photons in the system, 
the system cannot be classically simulated via truncation 
unless the level of truncation is very small. 




FIG. 3: Colour online. Simulating a lossy system via truncat- 
ing a classical simulation to m photons and discarding higher 
photon number terms. 



